30th International Cosmic Ray Conference 




ICRC*/ 



Merida, Mexico 



o : 
o , 

(N 

m 

<n : 



43 . 

Or 

i ' 
o 



> 

On 

00 
cn 

ON 

o 
l> 
o 



X 



Results from the Blazar Monitoring Campaign at the Whipple 10m Gamma-ray 
Telescope 

David Steele 1 for the VERITAS Collaboration 9 , and M. T. Carini 2 , P. Charlot 3 , O. 
Kurtanidze 4 , A. Lahteenmaki 5 , T. Montaruli 6 , A. C. Sadun 7 , M. Villata 8 
1 Adler Planetarium & Astronomy Museum, 1300 S. Lake shore Dr., Chicago, IL 60605, USA 
2 Dept. of Physics and Astronomy, Western Kentucky University, 1906 College Heights Blvd #11077, Bowl- 
ing Green, KY 42101, USA 

3 Laboratoire d'Astrophysique de Bordeaux, UMR 5804 LAB, 2 Rue de I'Observatoire, BP 89, 33270 
Floirac, France 

4 Abastumani Astrophysical Observatory, 383792 Abastumani, Georgia 
Mets&hovi Radio Observatory, Metsahovintie 114, FIN-02540 Kylmala, Finland 
6 Dept. of Physics, Univ. of Wisconsin-Madison, 1 150 University Ave, Madison, WI 53706, USA 
7 University of Colorado at Denver, PO Box 173364, Denver, CO 80217, USA 
s Osservatorio Astronomico di Torino, via Osservatorio 20, 10025 Pino Torinese, Italy 
9 For full author list see G. Maier, "VERITAS: Status and Latest Results ", these proceedings 
contact: dsteele @ adlerplanetarium. org 

Abstract: In September 2005, the observing program of the Whipple 10 m gamma-ray telescope was 
redefined to be dedicated almost exclusively to AGN monitoring. Since then the five Northern Hemi- 
sphere blazars that had already been detected at Whipple are monitored routinely each night that they are 
visible. Thanks to the efforts of a large number of multiwavelength collaborators, the first year of this 
program has been very successful. We report here on the analysis of Markarian 421 observations taken 
from November, 2005 to May, 2006 in the gamma-ray, X-ray, optical and radio bands. 



Introduction 

Among active galactic nuclei (AGN), blazars are 
the most extreme and powerful sources known and 
are believed to have their jets more aligned with 
the line of sight than any other class of radio-loud 
AGN. Blazars exhibit a broad continuum extend- 
ing from the radio through TeV 7-rays. and are 
characterized by large, rapid, irregular amplitude 
variability in all accessible spectral bands. Their 
high level of variability makes long-term, well- 
sampled, multiwavelength observations of blazars 
very important for constraining and understanding 
their emission mechanisms and characteristic time- 
scales. 

In September 2005, the observing program of 
the Whipple 10 m imaging atmospheric Cherenkov 
telescope (IACT) was redefined to be dedicated 
almost exclusively to nightly blazar monitoring. 



Since that time the five Northern Hemisphere 
blazars that had already been detected at Whipple 
(Mrk 421, H1426+428, Mrk 501, 1ES 1959+650 
and 1ES 2344+5 14) have been monitored routinely 
each night that they are visible above an elevation 
of 60°. For the first time, this has provided the 
opportunity to obtain long-term and well-sampled 
VHE light curves of these highly variable objects. 
Part of the motivation for these observations was to 
provide a trigger for more sensitive VHE observa- 
tions of these AGN by the new generation of IACT 
telescopes (CANGAROO-III, HESS, MAGIC and 
VERITAS) and to provide baseline observations 
for similar observations with GLAST. 

Here we will focus on the results of the monitoring 
campaign on Mrk 421. At a redshift of z=0.031, 
Mrk 421 was the first extragalactic source of VHE 
7-rays to be discovered with ground-based tele- 
scopes [9]. Since then Mrk 421 has been stud- 
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ied extensively in the VHE regime and many pe- 
riods of intense variability have been observed [3, 
8, 1]. Mrk421 is a high-frequency-peaked BL Lac 
(HBL) object exhibiting a non-thermal, broadband 
spectral energy distribution (SED) with two peaks, 
one in X-rays and a second in the TeV regime. 
The spectrum has been successfully reproduced by 
models where the X-ray peak is due to electron 
synchrotron radiation, and the TeV emission arises 
as the same electron population cools by inverse- 
Compton up-scattering of either the synchrotron 
photons (Synchrotron Self-Compton; SSC) or am- 
bient photons (External Compton; EC) [7]. Like 
all such objects studied to date, the two peaks in 
the SED have been found to correlate with each 
other both in peak frequency and in strength. 

The call for full multiwavelength coverage in pre- 
vious campaigns was only invoked when large 
flares were observed in the TeV regime. For the 
campaign reported in this work, we make no "a 
priori" assumptions about flux levels so that we can 
look for correlations on longer times scales and un- 
der different source conditions. Here we present 
light curves of radio, optical, X-ray and 7-ray data 
and a search for correlated emission between spec- 
tral bands. 

Methods 

All the TeV 7-ray observations presented here 
were made with the 10 m Telescope at the Fred 
Lawrence Whipple Observatory (FLWO). During 
the period November 29, 2005 to May 29, 2006 
a total of 98.0 hours of data were taken on Mrk 
421. The data presented here were mostly taken 
in the Tracking mode and were analyzed using 
analysis procedures pioneered and developed by 
the Whipple Collaboration [4]. Cuts were derived 
which reject approximately 99.7% of the cosmic- 
ray background images while retaining over 50% 
of those generated by 7-ray showers. A large sam- 
ple of dark-field data (« 200 hours) spanning a 
similar zenith angle range was analyzed to estimate 
the background during AGN data runs. Although 
sensitive in the energy range from 200 GeV to 
10 TeV, the peak response energy of the telescope 
to a Crab-like spectrum during the observations 



reported upon here was approximately 400 GeV 
(subject to a 20% uncertainty). 

The X-ray data for this campaign were taken by the 
Rossi X-ray Timing Explorer (RXTE). Data from 
the Proportional Counter Array (PCA) and the All 
Sky Monitor (ASM) were used, but only the PCA 
data are presented here. The sensitive energy range 
of PCA is 3 - 15 keV, limited at the low-energy 
end by systematic effects in the detector, and at the 
high-energy end by low source counts above the 
background [6]. Mrk 421 was observed with the 
RXTE PCA from January 6 to May 31, 2006 in 
short daily exposures of 19 minutes, on average. 
Due to observing constraints, the the X-ray and 
7-ray observations were not always simultaneous, 
with the PCA observations starting on average 99 
minutes behind the Whipple observations. Analy- 
sis of the PCA observations followed the standard 
procedure outlined in the RXTE Cookbook. 

Eight optical observatories contributed data sets re- 
sulting in excellent optical coverage over almost 
the entire length of the campaign. They are the 
FLWO 1.2 m telescope; the Tenagra 0.8 m tele- 
scope in Tenagra, Arizona, USA; the Bradford 
Robotic Telescope in Tenerife, Canary Islands, 
Spain; the WIYN 0.9 m telescope at Kitt Peak, 
Arizona, USA; the 0.7 m telescope at Abastumani 
Observatory in Abastumani, Georgia; the 0.6 m 
Bell Observatory at Western Kentucky University, 
Bowling Green, Kentucky, USA; Bordeaux Obser- 
vatory in Floirac, France; and the 1.05 m REOSC 
telescope at Osservatorio Astronomico di Torino, 
Italy. 

The data from the various observatories were re- 
duced and the photometry performed indepen- 
dently by different analysts using different strate- 
gies. Most analysts used relative aperture photom- 
etry performed using standard routines in IRAF 
[10]. Magnitudes were calculated with respect to 
one or more standard stars of Villata [11]. Be- 
cause different observatories use different pho- 
tometric systems, and photometric apertures and 
the definition of the reported measurement er- 
rors for each nightly-averaged flux are inconsis- 
tent across datasets, we have adopted a simple ap- 
proach for the construction of the composite light 
curves whereby a unique flux offset is found for 
each spectral band of every instrument based on 
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Figure 1: Radio, optical, X-ray and 7-ray light curves of Mrk 421 in 2005-2006. 



nearly simultaneous observations between partici- 
pating observatories. 

The radio observations of this campaign come 
from three observatories taken at 10 frequencies: 
Metsahovi Radio Observatory, Kylmala, Finland 
(37 GHz), the University of Michigan Radio As- 
tronomy Observatory, Dexter, Michigan, USA 
(4.8, 8.0 & 14.5 GHz), and the Radio Astronomical 
Telescope of the Academy of Sciences (RATAN- 
600), Zelenchuksky, Russia (1, 2.3, 4.8, 8, 11 & 
22 GHz). Only the 37 GHz observations are re- 
ported here; the full set of radio data, (as well as 
a full X-ray and 7-ray spectral variability analy- 
sis) will be presented elsewhere. The observations 
from Metsahovi at 37 GHz have a bandwidth of 
1 GHz and were obtained using the 14 m-diameter 
radio telescope. 

The resulting light curves for the spectral bands de- 
scribed above are presented in Fig. 1. The source 
exhibits significant variability in all wavebands, 
with the level of variability generally higher in the 
energy regimes corresponding to the two peaks of 
the SED. 

Results & Discussion 

Of all the pairs of light curves investigated, the X- 
ray and TeV light curves demonstrate the strongest 
correlation, as demonstrated by plotting the TeV 
flux vs. X-ray flux as in Fig. 2. The correlation 
suggests that the particle populations responsible 
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Figure 2: TeV 7-ray flux vs. X-ray flux (3-15 keV). 

for at least part of the emission in the two energy 
regimes are collocated, as is the case in the SSC 
scenario. 

We searched for correlations between fluxes in 
other wavebands using the Discrete Correlation 
Function (DCF) [2]. The DCF gives the linear cor- 
relation coefficient (R) for two light curves as a 
function of the time lag (r) between them. Besides 
the X-ray and TeV light curves, only the DCF be- 
tween the optical and TeV light curves, presented 
in Fig. 3, warrants further investigation. The DCF 
indicates a possible correlation with the optical 
flux lagging the TeV flux by about 7 days. An 
elevated level of correlation is also seen with the 
optical leading the TeV by 25 to 60 days. 
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Figure 3: Discrete Correlation Function between 
the optical R-band and TeV band light curves. Pos- 
itive values of time lag correspond to the optical 
flux lagging the TeV flux. 

We investigated the significance of these possible 
correlations using simulated optical light curves 
with the same variability properties as the observed 
R-band light curve based on a Power Spectral Den- 
sity (PSD) derived from a fit to the first-order struc- 
ture function (SF). Assuming the process(es) re- 
sponsible for the optical variability are stationary 
over the period of observation, the slope of the 
PSD, (3, at a given variability frequency is related 
to the slope of the SF, a, by (3 = 1 + a [5]. 

After multiple realizations of the campaign using 
simulated, uncorrected optical light curves we as- 
sess the probability to have seen the observed cor- 
relation features due to chance. For each bin of 
each simulated DCF, we record the values of the 
correlation coefficient, R, its error, cfr, and S = 
R/(Jr. After many realizations, we build the dis- 
tribution of S seen in each bin of the DCF and 
compute the chance probability, pi, to have found 
S > S bs- To assess the likelihood of the observed 
DCF features as a whole, we construct two more 
distributions, £± = n«K where for £„ (£+) the 
product runs over all DCF bins with 60 d< r < 
(0 < r < 60 d). From these, we find the likelihood 
to have observed the optical precursor (optically- 
lagged) features by chance to be 20% (60%). Fur- 
ther investigation reveals that the lagged optical 
emission feature seen at r = 7 d is likely due 
to contamination from the R-band auto-correlation 
function, which also exhibits alternate peaks and 
valleys at multiples of 7 days. The possible reality 
of the optical precursor emission is more difficult 
to interpret since our method assesses all the corre- 
lation features on a given side of the DCF at once. 



Unfortunately, it would not have been possible to 
assess a particular feature "a posteriori" without an 
unknown trials penalty. 

Conclusions 

We have presented multiwavelength observations 
of the blazar Mrk421 taken during the first year of 
dedicated AGN monitoring with the Whipple 10 m 
Telescope. The TeV 7-ray and X-ray fluxes were 
well-correlated, consistent with previous observa- 
tions of objects of this class. There was an indi- 
cation that the optical flux may have led the TeV 
flux by 25 to 60 days, but a simulation reveals sim- 
ilar or stronger correlations arise by chance in one 
out of five simulated campaigns. We anticipate fu- 
ture long-term, multiwavelength observations will 
be important for further eludication of AGN emis- 
sion mechanisms. 
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